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It is well known that a symmetric soliton in coupled nonlinear Schrodinger (NLS) equations with 
the cubic nonlinearity loses its stability with the increase of its energy, featuring a transition into 
an asymmetric soliton via a subcritical bifurcation. A similar phenomenon was found in a dual-core 
system with quadratic nonlinearity, and in linearly coupled fiber Bragg gratings, with a difference 
that the symmetry-breaking bifurcation is supercritical in those cases. We aim to study transitions 
between symmetric and asymmetric solitons in dual-core systems with saturable nonlinearity. We 
demonstrate that a basic model of this type, viz., a pair of linearly coupled NLS equations with the 
cubic-quintic (CQ) nonlinearity, features a bifurcation loop: a symmetric soliton loses its stability via 
■ a supercritical bifurcation, which is followed, at a larger value of the energy, by a reverse bifurcation 

that restores the stability of the symmetric soliton. If the linear-coupling constant A is small enough, 
Q \ the second bifurcation is subcritical, and there is a broad interval of energies in which the system 

is bistable, with coexisting stable symmetric and asymmetric solitons. At larger A, the reverse 
bifurcation is supercritical, and at very large A the bifurcation loop disappears, the symmetric 
soliton being always stable. Collisions between solitons are studied too. Symmetric solitons always 
collide elastically, while collisions between asymmetric solitons turns them into breathers, that 
subsequently undergo dynamical symmetrization. In terms of optics, the model may be realized in 
both the temporal and spatial domains. 
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I. INTRODUCTION 



Twin-core nonlinear fibers and waveguides (alias couplers) have been a subject of considerable interest in nonlin- 
ear optics and, more generally, in nonlinear- wave theory since the pioneering works by Jensen and Maier 0, who 
• introduced the model of a dual-core nonlinear optical fiber. In particular, many papers analyzed dynamics of solitary 
waves in models of nonlinear couplers (see, e.g., review Q and Ref. Q), although such solitons have not yet been 
created in the experiment. 

It has been known for a long time 2] that a symmetric soliton in the model of the twin-core fiber, with its energy 
equally divided between the cores, becomes unstable when the energy exceeds a certain critical value. Then, it was 
found that this instability gives rise to a pitchfork bifurcation: a pair of new, stable asymmetric solitons (which are 
mirror images to each other) appear when the symmetric state loses its stability 6]. This bifurcation was studied in 
detail by means of numerical 0, Q and analytical (3 methods (the latter was based on the variational approximation, 
see Section 6 of review 0), with a conclusion that the bifurcation is slightly subcritical, according to the standard 
definition [l(| ■ Namely, new stable asymmetric states appear at a value of the soliton's energy which is slightly smaller 
than that at which the symmetric soliton becomes unstable, thus giving rise to the bistability (coexistence of stable 
symmetric and asymmetric solitons) in a narrow interval of energies before the pitchfork bifurcation. 

The studies of solitons in nonlinear dual-core systems were extended in various directions. First, the change of the 
character of the bifurcation was investigated in a practically interesting case of a dual-core nonlinear optical fiber 
with asymmetric cores |3.Hl|. Further, four-component solitons were considered in a more sophisticated model, which 
takes into regard two orthogonal polarizations of light in the dual-core fiber Then, bifurcations breaking the 

symmetry of two-component solitons were studied in a symmetric dual-core nonlinear fiber with the Bragg grating 
carried by each cor e Il3l . and in a system of parallel-coupled waveguides with the quadratic, (rather than cubic, 
X^)i nonlinearity 14] . In the two latter cases, the symmetry- breaking bifurcation was found to be supercritical, 
rather then subcritical (i.e., the destabilization of the symmetric soliton immediately gives rise to a pair of stable 
asymmetric ones, without the bistability). In addition, three-component solitons and their bifurcations were studied 
in a system of three linearly-coupled fibers forming a triangular configuration |l5l | , as well as in a similar configuration 
formed by three Bragg gratings 

A universal feature of all the above-mentioned systems is that the symmetric soliton is destabilized with the increase 
of its energy, as the self-focusing nonlinearity favors states in which a larger part of the energy is concentrated in one 
core (which is explained by the fact that the part of the system's Hamiltonian accounting for the nonlinearity tends 
to be lower in such an asymmetric configuration than in its symmetric counterpart). Past the symmetry-breaking 
bifurcation point, the asymmetry of the energy distribution between the cores increases monotonously with the further 
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increase of the total energy. 

However, the self-focusing nonlinearity may feature saturation with the increase of the power in some optical media. 
In particular, it was inferred from experimental data that nonlinearity of chalcogenide glasses [17| and some organic 
transparent materials is adequately approximated by adding a self-defocusing x^ (quintic) term to the self- 
focusing %( 3 ) one. The respective cubic- quintic (CQ) nonlinearity furnishes a paradigmatic example of the saturation, 
as well as of the competition between self- focusing and self-defocusing nonlinearities |Tflj . 

A model of a dual-core optical hber or waveguide with saturable nonlinearity is an interesting object in its own right, 
and may also be relevant for applications to all-optical switching. In particular, the symmetry-breaking bifurcation 
between continuous-wave (CW) states (rather than solitons) in a model of a coupler with saturable nonlinearity was 
considered long ago in Ref. |2(j, with a conclusion that the bifurcation is subcritical in this case. However, CW states 
in couplers with a self- focusing nonlinearity are always subject to modulational instability |2l| , therefore it is necessary 
to consider solitons, rather than the CW. Qualitatively new features may be expected in such systems. Indeed, since 
the spontaneous symmetry breaking of soliton states, observed at energy exceeding the critical value, is accounted 
for by the trend to concentrate a larger part of the energy in one core with the growing strength of the self-focusing, 
the saturation, which attenuates the self-focusing with the further increase of the energy, should eventually give rise 
to a reverse trend. Thus, a reverse bifurcation may occur in the system, restabilizing the symmetric soliton and 
eliminating asymmetric states. 

This work aims is to produce a full picture of bifurcations in the symmetric dual-core system with the CQ non- 
linearity. In Section II, we introduce the model, whose single free parameter is the constant A of the linear coupling 
between the cores. Main results are presented in Section III: if A is smaller than a certain value A max , the system 
features a bifurcation loop, first destabilizing and then restabilizing the symmetric soliton, as may be expected (see 
above). The direct (symmetry-breaking) bifurcation is always supercritical, while the reverse one is subcritical (fea- 
turing a large bistability region, which may be of considerable interest to applications - Q), except for a narrow 
interval near A = A max , where the reverse bifurcation is supercritical, without the bistability. The loop disappears 
at A = A max ; in the region of A > A max , the symmetric soliton is always stable, while asymmetric ones do not exist. 
Direct simulations demonstrate that the stability of all soliton branches exactly complies with what may be expected 
from general principles of the bifurcation theory [lOj . 

In Section IV, we additionally investigate collisions between moving solitons, which a natural issue to address as 
the Galilean invariance of the model allows one to generate moving solitons. We infer from simulations that collisions 
between stable symmetric solitons are completely elastic, while asymmetric solitons collide inelastically and separate 
after the collision in the form of breathers. Conclusions are formulated in Section V. 



II. THE MODEL 

Linearly coupled CQ nonlinear Schrodinger (NLS) equations for the amplitudes u and v of the electromagnetic 
fields propagating along the coordinate z in the two cores of an optical waveguide can be cast in a normalized form, 

iu z + u TT + 2\u\ 2 u — |w| 4 ?i + Xv = 0, (1) 
iv z +v TT + 2\v\ 2 v- \v\ 4 v + Xu = 0, (2) 

where r = t — z/Vq is the local time (Vq is the group velocity of the carrier wave). The model assumes anomalous 
group- velocity dispersion, GVD (otherwise it will generate no bright solitons), whose coefficient is normalized to be 
1. The coefficients in front of the self-focusing x^ an d self-defocusing x terms can always be normalized as set 
in Eqs. Q and the only remaining irreducible parameter of the model being the linear-coupling constant A (we 
define it to be positive), that determines a characteristic coupling length, z COU pi ~ 1/A (in real couplers, z C oupi takes 
values from ~ 1 mm up to ~ 1 cm, in physical units). In nontrivial soliton states (those which are essentially affected 
by the linear coupling), the characteristic nonlinearity and dispersion lengths, z n oniin and z disp, must match z coup \. 
The successful creation of solitons in fiber Bragg gratings, with the nonlinearity length < 1 cm [52, suggests that 
Znoniin for high-power pulses in silica fibers can indeed be brought down to the necessary size, in the range of 1 mm - 
1 cm (the nonlinearity length is determined by the core itself, rather than by its coupling to another core or grating 
written on its surface). On the other hand, while the intrinsic GVD of the fiber waveguide may not be sufficient to 
provide for the necessarily small Zdi sp , the effective GVD may be enhanced by several orders of magnitude by means 
of the same grating [22j| . 

In addition to the above-mentioned interpretation in the temporal domain, the model based on Eqs. an d @ 
may also be realized in terms of the light transmission in the spatial domain. In that case, the temporal variable 
r is replaced by the transverse coordinate x, and the equations govern the spatial evolution of amplitudes of the 
electromagnetic field in a pair of parallel tunnel-coupled planar waveguides, the terms u xx and v xx accounting for the 
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transverse diffraction of light (in the usual paraxial approximation) . Launching the beams with the transverse width 
Aa; ~ 20 wavelengths (if the latter is taken, as usual, in a ballpark of A ~ 1 fim; note that the paraxial approximation 
is relevant for Ax ~ 20A) will provide for a value of the diffraction length, Zdiffr = 4-7T (Ax) 2 /A, in the above-mentioned 
range of 1 mm - 1 cm, where it can be matched to the coupling length (for the parallel planar waveguides, z coup i is 
essentially the same as for the dual-fiber system). 

The starting point of analysis is a well-known exact soliton solution of the single CQ NLS equation [23j, to which 
Eqs. and reduce in the symmetric case: 

u = v = e ikz U aymm (T), 



2 (fc - A) 

Symm(T) ~ Vl + Vl-4(fc-A)/3cosh(2Vfc-Ar)' ( } 

where the propagation constant k takes values in the interval A < k < | + A. In more sophisticated settings, solutions 
for a single equation with the CQ nonlinearity can be found in a numerical form. An example relevant to the spatial- 
domain model is a solution for the semi-nonlinear equation, with an overall nonlinear coefficient (one in front of both 
the cubic and quintic terms) present in the region of x > and vanishing at x < 

Following the analysis of coupled NLS equations with the nonlinearity we look for asymmetric stationary 

solutions to Eqs. and J2J in the form of 

{u(z,t),v(z,t)} = e* kz {U(t),V(t)} . (4) 

If the functions U and V are complex, 

{U(t),V(t)} = {a(r)e^\b(T)e^} , (5) 

with real amplitudes and phases {a, b} and {<f>, ip}, a straightforward corollary of the equations is the following relation, 
valid for localized solutions (see details in Appendix): 

a 2 </> + b 2 ijj = 0, (6) 

with the overdot standing for d/dr. For symmetric solitons, with a = b and <f) = tp, Eq. JHJ implies, as usual, 
that all solutions are real. While the same is not obvious for asymmetric solitons, we demonstrate in Appendix 
that asymmetric complex solutions for U(t) and V(r) with a nontrivial phase structure cannot be generated by a 
bifurcation from the symmetric soliton J3J. For this reason, in what follows below we consider real stationary solutions 
for U and V, hence the substitution of expressions J3J in Eqs. QJ and (J2J leads to a system 

42 ' kU + XV + 2U 3 - U b = 0, (7) 
dr z 

d 2 V 

—5 — kV + XU + 2V 3 — V 5 = 0. (8) 
dr £ 

We solved Eqs. J7J) and iJSJ numerically by means of the finite-difference method. In this paper, we only aim to 
search for fundamental solitons, which correspond to even single- humped solutions to these equations. Stability of the 
analytical symmetric solutions © and asymmetric ones found in the numerical form was tested by direct simulations 
of Eqs. ||TJ and performed by means of the split-step method. 

Note that Eqs. and @ conserve the total energy (norm) of the solution, 

1 r +oc 

N =2] 00 [HT)\ 2 + \v(T)\ 2 ]dT (9) 

(in the spatial-domain optical model, the norm has the physical meaning of the total power of the beams), together 
with the momentum and Hamiltonian (that will not be used below). For exact solutions @, the norm is 



\/3 / V3 + 2y r k~^X ^ 
~2~ n \VS~2Vk = ^t 



N symm = -V In ( v v ) . (10) 



It is worthy to note that expression ljTU|l satisfies the condition dN/dk > 0, which, according to the known Vakhitov- 
Kolokolov (VK) criterion [25|], is necessary for stability of the soliton family (the validity of this criterion was es- 
tablished, in a general form, for single-component NLS equations [2^|). The VK criterion guarantees the absence of 
unstable modes of small perturbations with a real instability growth rate, but it does not tell anything about modes 
corresponding to complex eigenvalues. 
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III. THE BIFURCATION DIAGRAM 



Symmetry-breaking bifurcation points of Eqs. J7J) and |JSJ (critical values of k for given A) may be looked for by 
adding an infinitesimal antisymmetric perturbation, W(t), to the symmetric soliton, {U(t), V(t)} — U symm (T)±W (t) . 
The substitution of this in Eqs. Q and (jSJ and linearization with respect to W lead to a linear Schrodinger equation 
for W, 

-(X + k)W = -^- + U aB (r)W, (11) 

with an effective potential 

U cS (t) = C/ S 2 ymm (r) [5C/ 2 symm (r) - 6] . (12) 

The first problem is to find values of k at which an even localized state appears among solutions of Eq. 1|11JI . In 
the model with the cubic nonlinearity only, Eq. i|12|) yields the effective potential — 6(fc — A)sech 2 (Vfc — At), which 
is integrable, as is known from quantum mechanics. Therefore, the (single) bifurcation point in the model of the 
coupler with the cubic nonlinearity could be found exactly 2]. In terms of the norm of the symmetric soliton, it is 
^cubic = 2-y/2A/3. However, the effective potential in the CQ model, given by Eqs. (jl2(l and 0, is far from being 
integrable, therefore we will identify bifurcation points by numerical methods. 

Numerical solution of stationary equations JTJ and JSJ produces a sequence of bifurcation diagrams displayed in 
Figs. ^ and |21 In these diagrams, the soliton's asymmetry parameter, 

U 2 - V 2 

^ max y max o\ 

6 - TJ2 +V 2 > ^) 
w max 1 v max 

where C/ max and V max are the amplitudes of the two components of the soliton, is shown versus the norm N. Note 
that the soliton's propagation constant k varies along with N (not shown in the diagrams). 

As expected, Figs. Q] and [21 demonstrate bifurcation loops. The loop exists at A < A max w 0.44 (in fact, the full 
loop was found for A > 0.10, as for smaller A it was difficult to close the loop, due to its very large size; indeed, for 
A — > 0, the loop extends to N — > oo, as one then has an obvious stable asymmetric solution, with u ^ and v = 0, 
for any N). The direct (symmetry-breaking) bifurcation is observed to be always supercritical, while the reverse one, 
which closes the loop, is subcritical (giving rise to the bistability and concave shape of the loop, on its right side) up 
to A ss 0.40. In a narrow interval of 0.40 < A < 0.44, the reverse bifurcation is supercritical, and the (small) loop has 
a convex form. 

It is relevant to mention that a somewhat similar sequence of bifurcation loops was reported for (modulationally 
unstable) CW states, rather than solitons, in a symmetric dual-core model by Snyder et al. [2(| [see Figs. 10(a) 
and 10(b) in that paper]. However, that sequence was found in a model with a saturable, rather than cubic-quintic, 
nonlinearity in each core. 

The global picture of the bifurcations is additionally illustrated by Figs. |3and0] They display, versus A, the value 
of the asymmetry parameter <|13[l for the most asymmetric soliton generated by the bifurcation for given A, and values 
of norm N of the symmetric soliton at points of the direct and reverse bifurcations. 

Stability and instability of different branches of soliton solutions can be anticipated on the basis of general principles 
of the bifurcation theory [lOj : the symmetric solution becomes unstable after the direct supercritical bifurcation, and 
asymmetric solutions emerge as stable ones at this point; after the reverse bifurcation, the symmetric solution is stable 
again. In the case when the reverse bifurcation is subcritical and, accordingly, the bifurcation loop is concave on its 
right side, two branches of asymmetric solutions (one stable and one unstable) meet and terminate at each turning 
point through the saddle-node (alias tangent) bifurcation; accordingly, the branches of the asymmetric solitons which 
terminate at the reverse-bifurcation point are unstable. 

These expectations are fully borne out by direct numerical simulations. In particular, the simulations clearly 
demonstrate that a small symmetry-breaking perturbation triggers spontaneous rearrangement of presumably unstable 
symmetric solitons into stable asymmetric ones. The transition (not shown here) is accompanied by relatively small 
radiative loss. 

Simulations of the evolution of perturbed unstable solitons, corresponding to the intermediate branches in Figs. ^ 
and|21(i.e., in the case of the bistability), demonstrate that the unstable soliton, which has a choice to transform itself 
into either a still more asymmetric soliton, or a symmetric one, both of which are stable, clearly follows the latter 
route, even if the initial perturbation acted in the opposite direction, trying to make the soliton less asymmetric (not 
shown here in detail). In this case, the transition entails a larger radiative loss. 

The instabilities of symmetric and asymmetric solitons are qualitatively different: as said above, the symmetric 
solutions (Q are stable according to the VK criterion, i.e., their instability may only be oscillatory [accounted for by 
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FIG. 1: A set of bifurcation diagrams for symmetric and asymmetric solitons, in the plane (N,e), as found from numerical 
solution of Eqs. J7J and (|HJ at different values of the linear-coupling constant A. Stable and unstable branches of the solutions 
are shown by solid and dashed curves, respectively, and bold dots indicate bifurcation points. 
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FIG. 2: Continuation of Fig. □ 
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FIG. 3: The maximum value of the asymmetry parameter 1)13^ of solitons generated by the bifurcation, versus the linear- 
coupling constant A. For very small values of A, e max is not shown because of difficulties with convergence of the numerical 
scheme. 
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FIG. 4: Values of the norm (energy) of the symmetric soliton at which the direct and reverse bifurcations occur. The two 
curves merge and terminate at A = A max ~ 0.44. 



complex eigenvalue(s)], which is indeed observed in the simulations (the validity of the VK criterion was not proved 
for the present model, but it was rigorously derived for another system, including two linearly coupled equations with 
the cubic nonlinearity of opposite signs |27|). On the other hand, inspection of dependences N(k) for the asymmetric 
solutions demonstrates that, in all cases when these solitons are unstable, they have dN/dk < 0, i.e., they are expected 
to be VK-unstable solutions. The latter implies that the instability, being accounted for by a real eigenvalue, must 
grow without oscillations. The simulations support this expectation. 

Stability of all solution branches which are shown as stable ones in Figs. ^ and [21 has also been verified in direct 
simulations. In all these cases, it was found that even relatively large perturbations do not destroy the solitons. 
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FIG. 5: A typical example of the elastic collision between two identical in-phase symmetric solitons generated by means of 
expressions 114H with x — ilj f° r A = 0.32 and N = 1.106 (which correspond to k = 0.56). Only the u-component is shown 
here, as the picture in the M-component is identical. 



with an arbitrary real constant x generates a boosted solution differing from the original one by r — > r + 2\z , i.e., 
it moves at the (inverse) velocity c = — 2\ relative to the original state. The invariance suggests to study collisions 
between moving solitons (collisions between symmetric and asymmetric solitons in the model of a coupler with the 
cubic nonlinearity were investigated in Ref. |28j). Simulations of the collisions make it possible to draw the following 
general conclusions. 

First, collisions between stable symmetric solitons (which are obviously tantamount to collisions between solitons 
in the single-component NLS equations with the CQ nonlinearity) appear to be completely elastic, as shown in Fig. 
El The elastic collisions give rise to small shifts of the solitons' centers, in the same fashion as in collisions between 
solitons in the integrable cubic NLS equation (each soliton is additionally shifted in the direction of its motion). 

Second, collisions between stable asymmetric solitons are essentially inelastic. In this case, the solitons emerge from 
the collision with smaller velocity c, and in an excited state, i.e., as breathers, rather than stationary solitons. For 
example, the (inverse) velocities of both solitons drop from q n = ±0.4 before the collision to c out = ±0.2 after the 
collision, if the moving solitons were generated by means of expression Ijl4(l with \ = ±0.2, for A = 0.3 and N — 1.176 
(which correspond to k = 0.6). 

The breathers which appear after the inelastic collision are not persistent ones. As Fig. shows, the amplitude of 
the intrinsic oscillations of the breather gradually decreases, while the oscillation period increases. This may imply 
that, in the course of subsequent extended evolution, the breathers will eventually relax to stationary solitons (due to 
losses induced by emission of radiation), but detailed study of this stage of the evolution is beyond the scope of this 
work. It is noteworthy that, while the solitons generally keep their identity after the inelastic collision (the solitons 
which originally had a larger u- or v-component keep this difference after the collision), the post-collision oscillatory 
evolution clearly tends to symmetrize both pulses, so they may eventually settle down to a pare of identical symmetric 
solitons. 



IV. COLLISIONS BETWEEN STABLE SOLITONS 
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FIG. 6: The evolution of one of the breathers produced by the inelastic collision of two solitons in the above-mentioned example, 
with A = 0.3, N — 1.176, k — 0.6, and shove parameters x = ±0.2 (the second breather is a mirror image of the one shown 
here). The solid and dashed curves show the amplitudes of the u- and v-components of the breather. 

V. CONCLUSION 

In this work, we have extended the model of the nonlinear dual-core coupler by adding defocusing quintic nonlinearity 
to the usual cubic self-focusing term. In terms of optics, the model may be interpreted in temporal and spatial domains 
alike (i.e., in fiber couplers and in a pair of tunnel-coupled parallel planar waveguides). As a result, at not too large 
values of the linear-coupling parameter A, we observe a bifurcation loop for solitons: the supercritical bifurcation, 
which destabilizes the symmetric soliton and creates a pair of asymmetric ones, is followed by a reverse subcritical 
bifurcation that restores the stability of the symmetric soliton. In this case, the loop has the concave shape on its 
right side, and includes a large bistability region, which may be of interest for applications to all-optical switching. At 
larger values of A, the loop's shape becomes convex, and the bistability disappears. At still larger A, the symmetric 
solitons always remain stable, and asymmetric ones never appear. 

The stability of all the branches of the soliton states, verified by direct simulations, exactly complies with what may 
be predicted, on the basis of the shape of the bifurcation diagrams, by general principles of the bifurcation theory. 
Additionally, in those cases when the asymmetric solitons are unstable, their instability is correctly predicted by the 
Vakhitov-Kolokolov criterion (while the instability of symmetric solitons has a different nature, as they are unstable 
against oscillatory perturbations). The simulations also demonstrate that unstable symmetric solitons rearrange 
themselves into stable asymmetric ones, with small radiation loss. Unstable asymmetric solitons rearranged into their 
still more asymmetric stable counterparts. 

The Galilean invariance of the model makes it possible to study collisions between moving solitons. Symmetric 
solitons collide in a completely elastic fashion, while collisions between asymmetric solitons are inelastic: they lose, 
roughly, half of the velocity, and emerge from the collisions as breathers, that subsequently tend to symmetrize 
themselves and, possibly, settle down to stationary symmetric solitons. 

It is plausible that basic results reported in this work are common to solitons in double-core systems with various 
forms of saturable nonlinearity. It may also be quite interesting to consider similar problems for two-dimensional 
solitons, that can be realized as spatiotemporal "light bullets" in parallel-coupled planar waveguides with the 
intrinsic nonlinearity of the cubic-quintic type. 

Appendix: Absence of asymmetric complex solitons generated by a bifurcation from the symmetric soliton 

Here we aim to consider in more detail a possibility of the existence of stationary asymmetric soliton solutions in 
the complex form, as per Eqs. ®. Substituting the latter in the stationary version of Eqs. Q and @, one can easily 
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derive the following pair of equations, 



d 



= Xab sin ( 



tp) , — (b 2 ^ — — Aa&sin (0 — tp) . 



(15) 



and another pair of equations, that we do not need here in an explicit form. An obvious consequence of Eqs. i|15|) f° r 
localized solutions, with o(|r| = oo) = 6(|r| = oo) = 0, is relation @). 

If there is a bifurcation from the real symmetric soliton (J3J) that gives rise to complex asymmetric solitons with a 
nontrivial phase structure, then, in an infinitesimal proximity to the bifurcation point, where both the phase functions 
4>(t) and , 0( T ) an d deviations of a(r) and b(r) from expression © are infinitely small, Eq. © yields, in the lowest- 
order approximation, <f>-\-ip = 0, or, in other words, VK 1 ") = ~4>( T ) (it is trivial to remove a phase constant). Further, 
in the same approximation, an equation for an infinitesimal function </>(t) following from Eqs. I|15(l is 



^ + 2A fl4 (T) 



0. 



where the expression (0 is to be substituted for a, and 



T : 



dr' 



a 2 (r') 2 (Ar — A) 




- - sinh (2\/k - Xrj 



(16) 



(17) 



(T is a monotonous function of r, varying together with it from — oo to +00). Equation (|16|l with A > is a linear 
Schrodinger equation with a coordinate T and potential barrier (rather than a potential well), which, obviously, 
gives rise to no localized eigenfunctions, hence it cannot generate any bifurcation that would be signalled by the 
appearance of such a solution (if one sets A < 0, the conclusion will be the same, as in this case one will start with 
the antisymmetric soliton). 
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